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ABSTRACT 



The Galactic bulge source MOA-2010-BLG-523S exhibited short-term devi- 
ations from a standard microlensing lightcurve near the peak of an Aj^ax ~ 265 
high- magnification microlensing event. The deviations originally seemed con- 
sistent with expectations for a planetary companion to the principal lens. We 
combine long-term photometric monitoring with a previously published high- 
resolution spectrum taken near peak to demonstrate that this is an RS CVn 
variable, so that planetary microlensing is not required to explain the lightcurve 
deviations. This is the first spectroscopically confirmed RS CVn star discovered 
in the Galactic bulge. 

Subject headings: stars:spots - stars:variables:other - gravitational lensing - plan- 
etary systems 
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1. Introduction 



High-magnification microlensing events provide a powerful tool for planet detection, 
partly because planets are more likely to perturb these events and partly because their high 
magnification (hence high signal-to-noise ratio) allows even very small perturbations to be 
detected. However, non-microlensing flux variations are also enhanced in these events. In 
this paper we report on the discovery of an apparent planet candidate that turned out instead 
to be a highly magnified active star and discuss methods by which we identified and excluded 
this interloper. 

Stars are intrinsically variable, and star spots can induce substantial light curve varia- 
tions in cool stars. However, for most G and K dwarfs this variability is manifested at a low 
level because magnetic activity decays quickly with age. Late M dwarfs can remain active 
for a Hubble time, but they are faint and will not be common microlen sing source s. There 
is, however, an important sub-population of highly active RS CVn stars fjHalllll976[ ) that are 
intrinsically luminous. 



Magnetic activity is governed by the Rossby number (INoyes et al.lll984l ). Rq = P/tc, 
where P is the rotation period and tc is the convective overturn timescale. Greater rotation 
(smaller P) induces faster buildup of magnetic fields. Deeper convection (bigger tc) permits 
the fields to build up for a longer time before they propagate to the surface. In the Rq 
regime of interest here, the ob served rms photometric variability Ar is a very steep function 
of 1/Ro JHartman et al.lboOQh 



OC Rq ] 



Rq = —■ 



As stars leave the main sequence they will develop deep surface convection zones as they 
become cooler, but they will also be expand substantially and slow down due to angular 
momentum conservation. Hence, some special circumstance is required to induce or permit 
relatively rapid rotation. There are three potential mechanisms. First, a K dwarf may find 
itself in a close binary (either by birth, or through 3-body interactions), and thus be spun up 
by tides. Second, an F or G dwarf may find itself in a wider binary that is not initially tidally 
interacting. But as the dwarf evolves into a K subgiant, its expanding radius enables tidal 
interactions with its companion that then spin up the subgiant. Finally, stars in a narrow 
range of masses, 1.25 Mq ^ M < 1.5 Mq (typically F dwarfs) can spend most of their lives 
spinning fairly rapidly because of their shallow convection zones and so are still spinning 



^Based on observations made with the European Southern Observatory telescopes, Program ID 85. B- 
0399(1) 
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when they evolve into K subgiants. The lower mass limit is required for fast rotation to 
survive. Above the upper mass limit, stars evolve so rapidly through the Hertzsprung gap 
that they spend almost no time as subgiants. The resulting single-star RS CVn subgiants 
therefore span a narrow range of ages, 7 Gyr ^t>3 Gyr. 

Here we report the detection of the first spectroscopically confirmed RS CVn star in the 
Galactic bulge. The detection was beyond serendipitous. It resulted from intensive spectro- 
scopic and photometric observations of an extremely rare high-magnification microlensing 
event of a subgiant source. Only about 1 bulge subgiant per 100 million is so magnified 
each year. The intensive photometry was carried out to find planets (orbiting the lenses), 
while the high-resolution spectrum of MOA-2010-BLG-523S was obtained to study chemical 
abundances of bulge dwarfs and subgiants. 

MOA-2010 -BLG-523S is a subgiant, with a temperature T ~ 5123 K and surface grav- 
ity log^ = 3.6 JBensbv et al.l l201ll. l2013llR. As such, either the second or third mechanism 
of forming RS CVn stars should apply. That is, it is either in a binary that was "tidally 
activated" by the growth of the primary as it evolved along the subgiant branch, or it is an 
isolated, retired F dwarf. The mere existence of an isolated RS CVn star would be evidence 
for intermediate- age bulge stars. Of course, with just one detection, one could not make a re- 
liable estimate of the fraction of bulge stars that are of intermediate age. But there are other 
lines of ev idence for such a p opulation, including age estimates of microl ensed dwarfs and 



subgiants (IBensbv et al 



2002; van Loon et al. 



2003 



20111) and asymptotic gian t branch (AGB) stars (ICole fc Weinberg 



Uttenthaler et al.ll2007l ). Thus, it would be of considerable inter- 



est to distinguish between the single-star and binary-star scenarios. Unfortunately, we find 
that both scenarios are plausible, given the available evidence, and so no definitive statement 
can be made regarding a putative intermediate-age population. 

A major focus of the present paper is the secure identification of the microlensed subgiant 
as an RS CVn star. However, the process of this discovery is of independent scientific interest. 
The event became a focus of attention because of deviations from standard microlensing 
seen over the peak. The J-band lightcurve was quite well fit by a planetary model, and 
hence was far "along the road" to being published as a microlensing planet, in which case 
it would have been only the 14th such planet. It was really only very small discrepancies 
that led to the gradual unraveling of this picture and the recognition that the deviations 



^The stellar parameters quoted in this work are taken from Bensby et al.l (|2013l ). These values are slightly 
revised from the ones originally given by iBensbv et al.l (|2011l ) . Because we discuss the history (Section |6]) 
of how MOA-2010-B LG-523S was recognized to be an RS CVn star, we report here, for completeness, the 
Bensby et al.l (|2011l ) parameters that were available at that time: T = 5250, \ogg = 4.0, [Fe/H]= +0.1, 
^ = 2.1kms-i. 
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at peak are most likely due to magnified star spots rather than a planet orbiting the lens 
star. The fact that irregular variability due to spots can be fit by planetary microlensing 
is sobering. As we discuss in Section El it implies that great care is required to securely 
identify microlensing planets in high-magnification events for cases of low-amplitude signals 
that lack clear microlensing signatures. 



Observational Data 



Microlensing event MOA-2010-BLG-523 (RA,Dec) = (17:57:08.9, -29:44:5 8) (lb) = 
(0.59 , -2.58) was alerte d by the Microlensing Observations in Astrophysics (MOA) (IBond et al 
200ll : ISumi et allboill ) collaboration at UT 08:46, 21 Aug 2010, and again 26.5 hours later 
as a potential high-mag event that would peak at Amax ~ 70 in four hours. In fact the event 
continued to rise for another 36 hours, which triggered much more intensive observations. 
At UT 16:51 23 Aug, the Microlensing Follow Up Network (/iFUN) issued a high-mag alert, 
predicting a peak at UT 02:00 to 04:00 and on this basis contacted the VLT bulge-dwarf 
spectroscopy group, advocating observations in that time interval. At the same time /iFUN 
organized its own continuous photometric observations using the 1.3m SMARTS telescope 
at CTIO to begin shortly after twilight. Very importantly in the present context, these 
observations were carried out with the ANDICAM camera, which is equipped with an opti- 
cal/infrared dichroic, so that it can take images simultaneously in, e.g., J-band and if-band. 

While the prediction of peak time turned out to be correct, VLT was unable to observe 
the event exactly when requested due to a conflict with technical activities, but did make a 
2 hour exposure (split in 4x30- min) with UVES on VLT beginning nea r twilight (UT 23:56 ). 
The main information on this spectrum has already been reported by lBensby et al.l (120111 ). 



There are two other v ery important data sets com i ng fr om the Optical Gravitational 
Lens Experiment (OGLE) (jUdalski et al.lll994l : lUdalskil |2003| ). The event itself was moni- 
tored by OGLE-IV, which began operations in March 2010. However, during 2010, OGLE- 
IV was in commissioning phase and so did not issue alerts. The data were first reduced in 
November 2010. Unfortunately, the target falls in a gap between chips in the new 32-chip 
OGLE-IV camera, meaning that the target was captured only when small pointing errors 
moved the target onto a chip, which occurred about 1/3 of the time. Due to the high quality 
of OGLE data (and despite the reduced coverage), it was already evident that the source 
was a low-amplitude variable and indeed it was checked (and confirmed) at the time of the 
image reductions that these variations were not due to chip-edge effects. Hence, AU had 
already suggested at this time that "the analysis of this object may be more complicated 
than expected". 
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The target also appears in OGLE-III, which took microlensing data from 2002-2009. 
And in addition, it was in a field that was the subject of a special high-cadence 46 day 
campaign in 2001 whose aim was to find transiting planets, during which it was observed 
786 times. 

In addition there were several other data sets, which in particular define the falling wing 
of the lightcurve extremely well. These include the RoboNet 2.0m Faulkes North Telescope 
(SDSS-z) in Hawaii, the PLANET 1.0m Canopus Telescope (J) in Tasmania (Australia), 
the PLANET 0.6m telescope (/)in Perth, Australia, and the following //FUN telescopes: 
Auckland 0.4m (/), Farm Cove 0.36m (unfiltered), Kumeu 0.36m (/), Molehill 0.3m (unfil- 
tered) (all in New Zealand), The 0.6m University of Canterbury B&C telescope intensively 
observed both wings of the light curve. Like the MOA 1.8m telescope, it is located at Mt. 
John, New Zealand. Finally the MiNDSTEp 1.5m telescope (/) in La Silla, Chile obtained 
data including a few points over the crucial peak region. 



3. Microlensing Analysis (Simple Version) 

As discussed in Section [6l a complete analysis of the microlensing event MOA-2010- 
BLG-523 is complicated by spots on the surface of the source (called MOA-2010-BLG- 
523S). However, it is possible to derive reasonably robust estimates of all the microlensing 
parameters required to constrain the source properties without detailed modeling of these 
complexities. 

We begin by simply excising the data within 0.8 days of the peak and fitting the rest of 
the lightcurve fiux F to the standard Einstein-Liebes-Refsdal-Paczyhski (1936, 1964, 1964, 
1986) 5-parameter form 

F{t) = f,AHt]) + h: u^=^l^}^±A, A=^^^ (2) 

MV M + 4 

Here A is the magnification, u is the projected source-lens separation in units of the Einstein 
radius, uq is the impact parameter, to is the time of closest approach, is the Einstein 
crossing time, /s is the source fiux, and /b is any blended fiux that does not participate in 
the event but is within the same point spread function (PSF) as the source. If there is more 
than one observatory, then each requires its own (/s, /b). We find 

to = 5432.603 ± 0.002; = 18.5 ± 0.5 day ,uo< 0.002 (3) 

and for the OGLE observatory 

/s = 19.33 ± 0.03; ^ = 0.03 ± 0.03. (4) 
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(All times are given in HJD' =HJD-2450000). 

Inspection of the relatively flat-peaked lightcurve shows that the lens crossed directly 
over the source and that the source crossing time is (crudely) of order t^, ~ 0.15 day, implying 
a source size (normalized to the Einstein radius) p = t*/tE ~ 0.008. Hence, because uq ^ p, 
(and noting that A — )■ u^^ for -u ^ 1) we can approximate the peak predicted magniflcation 

as 



A. 



1 + 



^-i)r 



(5) 



where (r") is the nth moment of the source surface brightness, and where we have assumed 
a linearly limb-darkened (and unspotted) source in making the evaluation, in which case the 
moments can generally be evaluated 

^ (3/2)!(l + n/2)! 



in 



n/2 + 1 



(3/2 + n/2)! 



(6) 



Here F is the "natural" form of the linear limb- darkening coefficien t, defined by surface 
brightness S{r) oc 1 - r[l - (3/2)(l - (r/p)2)i/2] f lAlbrow et al1ll999l ). It is related to the 
standard form u hj T = 2u/{3 — u). It is more "natural" in the sense that there is no net 
fiux associated with the limb-darkening term, which results in simpler formulae when written 
in terms of F. This includes not just the moment equations (j6]), but all formulae without 
exception. For examp le, the limb- darkening term in the standard formula for ellipsoidal 
variation (IMorrisI 119851 ) . (15 + u)/{3 — u) becomes simply (5 + 3F). 



We adopt Tj = 0.4 77 from I Clare tl ( 120001 ). by applying the stellar parameters measured 
bv iBensbv etaP fcoilh : T = 5123 K, [Fe/H]=+0.06, log^ = 3.6, ^ = 1.68 kms'^ Hence, 
^max = 2.17/ p. We evaluate Amax by taking the ratio of observed fiux at peak to the fit value 
of fs, and get very nearly the same answer, whether using the average of the two OGLE peak 
points, or a median estimate of CTIO near-peak points: Amax = 265. We thereby derive 



p = 111 I A, 



0.0082 ± 0.0003 



(7) 



where the error is derived from the 3% error in /g and a 3% error in the peak fiux due to 
spots. 



4. Observational Properties of MOA-2010-BLG-523S 



4.1. Baseline Variability 



As we will argue below, the rms variability of the source is about 3%. This is to be 
compared with the photometric errors, which are typically close to 10%. If the source were 
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a strictly periodic variable, then the period could easily be identified by folding the light 
curve, despite the low signal-to-noise ratio (S/N) of individual points. The situation is more 
complex for a quasi-periodic variable (as would be expected for a rotating spotted star). We 
are therefore quite fortunate that the source lay in a 2001 OGLE transit-campaign field, 
which was observed 786 times on 32 separate nights during a 46-day window. Binning the 
data by day, we therefore achieve errors of 0.02, which is comparable to the amplitude of 
the signal. The result is shown in Figure [1] The lightcurve gives the clear impression of 
variability with a period of order 12 days. 

We then use all the OGLE-III data to test for a quasi-periodic signal. If this is a 
spotted star, we expect that the underlying physical mechanism (rotation of the star) will 
be strictly periodic, but that the phase of the variations will drift over time as spots appear 
and disappear. As discussed above, except during the transit campaign, we are compelled 
to fold the data to pick up any signal at all. On the other hand, if we fold data over an 
interval that is too long, the result will suffer from destructive interference between different 
phase regimes. We therefore consider separate fits to the data for each of the nine seasons, 
2001-2009. In each trial, we hold the period fixed at a common value for all seasons. Hence, 
there are 28 parameters [Period -|- 9 x (phase, amplitude, zero-point)]. At P = 10.914 ±0.055 
days, there is an improvement of Ax^ = 69 relative to a fit for constant magnitude in each 
season (9 parameters), i.e., 19 fewer parameters. See Figure |2l 

We find that the phases are not consistent from one season to the next, suggesting that 
the variations are not strictly periodic. To further test this, we fit for a single phase and 
amplitude together with a zero-point offset for each season. This produces an improvement 
(relative to no periodic variations) of only 30 for 3 dof. Clearly the quasi-periodic variations 
are favored over strictly-periodic variations. 

4.2. Source is the Variable 

Faint sources in crowded fields are usually blends of several stars rather than discrete 
sources. And, of course, for microlensing events there is guaranteed to be at least one 
additional star along the line of sight in addition to the source, namely the lens. Hence, 
observing baseline variations does not in itself prove that the source is variable. However, 
from the microlens fit presented in Section [31 we know that the blend is at least 15 times 
fainter than the source. Thus, if it were responsible for the ~ 3% variations seen at baseline, 
it would itself have to vary at the > 50% level on ~ 11 day timescales. Such stars are 
extremely rare. Moreover, the chance is remote that one of these would happen to align 
with a source that (from other evidence we will present below) is expected to be variable. 



Therefore we conclude that it is MOA-2010-BLG-523S that is varying. 



4.3. Calcium H&K Emission 



Figure [3] shows the region of the calcium H&K lines in the UVES spectrum taken by 



J3; 
all 



Bensby et al.l ( 1201 ll ) near the peak of the event. The emission is extremely strong. We 
measure S'hk = 0.79 A by taking the ratio of the flux in these lines to the mean "continuum" 
i n the neighboring "V" and "R" regions (see Fig. [3]). For comparison llsaacson fc Fischer 
(I2OIOI ) found only 3 cases of comparable or larger S'hk among 234 "subgiants" in their 
survey of field stars. See their Figu res 11 and 12. We will dis cuss these in Section ESI below, 
but for the moment note that the llsaacson fc Fischerl (120101 ) stars are substantially redder 
and more luminous than MOA-2010-BLG-523S. 



4.4. Microturbulence Parameter ^ 

Figure H] shows the microturbulence par ameter £ plot t ed against temperature for 26 



microlensed dwarfs and subgiants as found by lBensby et al.l (l2010l . I2OIII ). MOA-2010-BLG 



523S has one of the largest ^. Moreover, it is well above the upper envelope of points on the 
low-temperature part of the diagram. This high "microturbulence" may reflect real turbulent 
motions on the surface of the star (as would be expected for an active star), but may in part 
reflect rotational motion. Since microturbulence represents a Gaussian velocity distribution 
that adds to line is quadrature with other effects, like instrumental resolution, unmodeled 
rotational motion will contribute to as 

. .2 {r^) (vsiniV 1 -0.21, ,2 , , 
Ae='^^^-^ = ^^{vsimf {u^p) (8) 

where vsmi is the projected rotational motion. 

In fact. Equation ([H]) applies to sources that are not differentially magnified, which is 
of course the usual case, but not the present one. If the lens were directly aligned with the 
source, then 

A ..9 (f^) (fsini)^ 1 — 0.3r, .0 , . 

= 2 — 6 — ^'^'"'''^ = 

i.e., roughly 2/3 of the non-differentially magnified case. For the actual geometry at the time 
of VLT spectra and J-band limb-darkening, we find below that A^^ = 0.2(f sini)^. Hence 
the measured ^ places an upper limit on vsini, 

vsmi< = 3.8 km s"^ (10) 
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4.5. Lithium 



In principle, it is possible to produce an isolated rapidly spinning subgiant (hence, an 
isolated RS CVn star) in an old population via stellar mergers. For example, a 10-Gyr solar 
mass star could begin evolving off the main sequence and swallow a smaller star, say 0.3 Mq, 
that had been its companion. This would both spin up the cannibal ar id pro v ide fr esh fuel 
to extend its life. The mass would be raised above the break in the iKraftI (jl970l ) curve, 
so that the star would not substantially spin down during its extended life. It would then 
evolve along the subgiant branch in a manner simila r to any other 1.3 M p) star. However, 
this scenario is ruled out in the present case because iBensby et al.l ( 120111 ) detected lithium 
with abundance log e (Li) = 1.6. Essentially all lithium would have been d estroyed if there 
had been a stellar collision (IHobbs &: Mathieulll99ll : lAndronov et al.ll2006l ). Thus, if MOA- 
2010-BLG-523S could be shown to lack companions, it would be of intermediate age. 



4.6. Radial Velocity 

The fraction of microlensing events toward the bulge whose source stars lie in the bulge 
(as opposed to the foreground disk) is > 95%. This is primarily because the optical depth to 
lensing is much higher due to the higher column of lenses. But this effect is also compounded 
by the fact that there are simply more bulge sources in these fields compared to disk stars. 
Nevertheless, if a source is weird in some way, its weirdness may be intrinsically connected 
with it being one of the small fraction of disk sources. This possibility is especially relevant 
in the present case because the disk is known to harbor a population of youngish subgiants, 
whereas the bulge is not. 



The source radial velocity (RV), Vr = +97.3kms~^ (IBensby et al.ll201ll ). makes it highly 
unlikely that it is in the disk because the expected value for disk stars is fr-,disk = +10 ± 
34kms^^ (compared to +10 ± lOOkms"^ for the bulge). 



4.7. Source Size 



Bensby et al.l ( 1201 ll ) derive an equivalent {V — I)o = 0.86 color from their spectroscopic 



solution (primarily from the temperature, but also taking account of the metallicity and 
gravity). We find from the microlens solution in Section 12] that the unmagnified source flux is 
A/ = 3.18 mag fainter than the clump. From the color-magnitude diagram of the neighboring 
field, there appears to be very little differential red dening. Hen c e A I ~ A/q. Based on 



the measured metallicity distribution of bulge stars, iNataf et al.l (|2012| ) estimate that the 
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absolute magnitude of the clump is M/^ci 
source is 

Mr s = Mr d + AJ - 



-0.12. Therefore, the absolute magnitude of the 



5 log 



D. 



3.06 - Slog 



cl 



where the last term i s the ratio of th e distances to the source and the clump. 



Bessell fc Brett 



(11) 

We then apply 



standard techniques (lYoo et al. I l2004h to evaluate the sou rce radius, first using 
(119881 ) to convert {V — I) ^ (V — K) and then using iKervella et al.l (120041 ) to obtain the 
/C-band surface brightness from the {V — K) color. Finally, we find 



IL = 2.15 R 



D, 



:i2) 



cl 



Note in particular that this derivation is independent of any assumption about the Galacto- 
centric distance Rq or the geometry of the Galactic bar, etc. Subgiants would be expected 
to have R^ > 2Rq. Hence, the source cannot lie substantially closer than the bulge because 
it would then be too small to be a subgiant (as indicated by its spectroscopic gravity). 



4.8. Consistency of Spectrum with Stellar Rotation 

If the source is rotating with a period P = 10.9 days, as seems indicated by the quasi- 
periodic variability seen in Figures [1] and |2l then the surface velocity is f = 2ttRs/ P = 
lO.Okms~^(-Rs/2.15-R0). The upper limit f sin i < 3.8 km s^^ (Eq. [10]) then implies i < 22°. 
This is a plausible value since randomly oriented stars will be uniformly distributed in cosi. 
That is, 7% of stars have i < 22°, which is small but not implausibly so. 



4.9. Consistency with Maoz-Gould Effect 



Maoz fc Gouldl (Il994j ) predicted that microlensing of rotating stars would generate an 
apparent RV shift, which would change during the course of the event. Their principal point 
was that the magnitude of this effect falls off only linearly with relative source-lens separation 
z = u/p, compared to the quadratic fall-off of photometric effects: 



Av 



(r^) V sin i 



1 - o.2r 



sm I 



-V sm I sm I 



i: 



(rO) 2z Az 
where is the angle between the source-lens separation and the projected spin axis. 



(13) 



The iBensby et al.l ( 120111 ) spectrum is actually composed of 4 30-minute exposures, cen- 
tered on HJD-2455432 = (0.510,0.531,0.552,0.573). Because the lens came very close to 
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source center, we will adopt z = {t — to)/t^. As we discuss in Section [6l to is not very accu- 
rately predicted by the lightcurve with the peak data removed and is actually approximately 
to = 5432.66 (compared to to = 5432.60 found in Section [3]). Therefore, at the four epochs, 
z = (1.00,0.86,0.72,0.58). We find numerically that the pre-factor in Equation f lT^ at these 
four epochs is (0.294,0.333,0.311,0.266). Thus, the maximum predicted relative shift is only 
0.067f sini sin0 < 0.31 km s~^. Based on cross correlation, the four spectra are consistent at 
this level. 



5. Possible Local Analogs 

We search for local analogs of MOA-2010-BLG-523S in order to better understand its 
nature and, to this end, begin by summarizing its characteristics. 



5.1. Summary of Characteristics 



From iBensby et al.l ( 120111 . |2013| ). we know the temperature, iron abundance, gravity, 
microturbulence, and lithium abundance: T = 512 3=b98K, [Fe/H]= .06 ± 0.07 logg = 
3.60 ± 0.23, ^ = 1.68 ± 0.20, e(Li) = 1.64 ± 0.10. iBensbv et aP J2OI1I ) also remark that 



the source has high sodium, which they note could be "fixed" by making it 500K hotter or 
increasing log q by 1 dex. However, they investigate these possibilities and reject them. While 
Bensby et al.l (120111 ) did not take account of differential magnification in their analysis, the 
im pact of such dif f erenti al magnification is quite small. For example, "Profile 35" considered 
by Ijohnson et al.l ( 20101 ) had much stronger differential magnification, but this affected the 
temperature by only 20K (see their Figs. 3 and 8). 



The baseline variability analyzed in Section 14.11 is best modeled as having constant 
period P = 10.914 days, but variable phase over 9 years, as would be predicted for a spotted 
star. We find variability amplitudes in these seasons (in mmags) of 20 ± 4, 40 ± 8, 25 ± 10, 
26 ± 9, 53 ± 9, 28 ± 7, 13 ± 10, 27 ± 9, 83 ± 19. Hence, a median of 0.027 mag, implying an 
rms variability of 2%. This is actually a lower limit, since each season's variability measure 
can be impacted by destructive interference between spot cycles at different phases. 

Finally, in Section 14.31 we measured calcium H & K emission of 5hk = 0.79 A. 
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5.2. Comparison by Rossby Number to M37 Sample 



From stellar models, we find that the convective overturn timescale for a T = 5123K 
subgiant is 3 times longer than for the Sun, while the measured period is 2.3 tirnes sho rter. 
Hence the Rossby number is 7 times higher. From Figure 17 of iHartman et al.l (|2009[ ). we 
observe that stars in M37 with similar Rq — 0.3 have rms variability in the range of 1-6%. 
Hence, the observed variability is quite consistent with locally observed stars. 



5.3. Comparison to Isaacson & Fischer Sample 



As discussed in Section 14. 3[ llsaacson fc Fischer! ( 120101 ) found only three stars with com- 
parable or greater Shk in their sample of 234 subgiants. These are Hipparcos stars HIP 
5227, 8281, and 97501, which have {V - K) = 2.43, 2.56, 2.59. They are thus consider- 
ably redder than M OA-2010-523S, which has (V" — /)o = 0.86 estimated from its spectrum 
( iBensby et al.ll201ll ). corresponding approximately to (V — K)n = 1 - 92. T hey are also about 
1-1.5 mag more luminous than the subgiants in the iBensby et al.l ( 1201 ll ) sample shown in 



Fig ure IH Indeed stars of this color and luminosity would not be deliberately selected for 
the IBensby et al.l (120111 ) "dwarf and subgiant" program, and are excluded from the analysis 
if they are observed by accident (Johnso n et al., in prep ) . Bas ed on their V/K photometry, 
Hipparcos distances, combined with the iKervella et al.l (120041 ) surface brightness relations, 
these three stars have radii of 5.4, 5.9, and 7.0 Rq, respectively. 

All three of these stars are spectroscopic binaries and at least the first two are broadly 
consistent with being t idally synchronized. Their binary periods are respectively 27.3 and 
30.1, days (Eker 20081). which would imply surface velocities of lOkms"^ in both cases, 
while llsaacson fc Fischerl ( l2010l ) report v sini of 14 and 6.2 kms"'^, respectively. However, D. 
Fischer (2011, private communication) notes that the profile of the first star is contaminated 
by lines from a companion, which she estimates to broaden the v sin i determination by 20- 
30%, making both stars quite consistent with tidal synchronization. D. Fischer also notes 
that the third star (Hip 97501) is a clear double-lined spectroscop ic binary, so that its binary 
nature is not in doubt even though the 3 RV measurements by llsaacson fc Fischerl (|2010[ ) 
show a scatter of only 0.2kms~^. 

A plausible scenario for these stars is that their moderately close companions only 
started to spin them up when they began to expand their envelopes as they approached the 
giant branch. In particular, for the first two, their known periods (P ~ 30 days) are too 
long for tidal synchronization while the stars were on the main sequence. However, as they 
evolved along the subgiant branch, they were clearly tidally spun up, as evidenced both by 
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their v sin i and their calcium H & K activity. 

Because MOA-2010-BLG-523S is a much smaller star, a much closer companion would be 
required to tidally couple with it. On the other hand, its period is much shorter. Since tidal 
amplitudes scale oc R^/P^, a companion in an 11 day period could provide tidal interactions 
only a factor ~ 2 smaller than these local analogs. Because MOA-2010-BLG-523S's period 
is shorter and its phase of subgiant evolution is longer, it would have many more periods 
to tidally synchronize. Hence, tidal spin-up by a binary companion is a very plausible 
explanation for its variability and the strength of its calcium lines. 



5.4. Comparison to Kepler Sample 



Chaplin et al.l ( 1201 ll ) find a dramatic drop in the interval 5150K< T <5400K in the 
fraction of Kepler asteroseismology targets for which they can measure oscillations. These 
show variability in the range 0.1-10 mmag, which is modestly higher than neighboring tem- 
perature ranges. See their Figure 1. The most plausible interpretation is that subgiants in 
this temperature range preferentially acquire spots that physically interfere with the propa- 
gation of stellar oscillations. (Oscillations in dwarf stars at these temperature ranges would 
be undetectable in any case.) This cannot be due to close binary companions because the 
phenomenon is nearly universal, whereas only a few percent of stars have such close com- 
panions. Rather, the physical mechanism must be that as single stars evolve redward on the 
giant branch, their convection zones deepen, so tc increases, and they become more spotted. 
After they pass through the most-affected temperature range, they expand rapidly, thus in- 
creasing their moment of inertia and so slowing their rotation. Of course, the more rapidly 
they are rotating at the outset, the higher the level of activity, but the increase in activity 
at this temperature range is nearly universal. 

Because the temperature T = 5123 K of MOA-2010-BLG-523S is at the edge of this 
affected range, it is also a plausible candidate to be a non-binary active subgiant. 



6. Birth and Death of a Microlens "Planet" 

Due to its predicted high-magnification, MOA-2010-BLG-523 was monitored almost 
continuously over peak by the 1.3m SMARTS telescope, although there were short gaps to 
check on another, possibly interesting, event. These data, by themselves, display a significant 
"bump" near peak. Moreover, the time of the observed peak is asymmetrically offset from 
that expected based on the MOA data (roughly a half day on either side of peak) by about 
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1.5 hours. Such bumps and asymmetries are just the type of features we look for to identify 
planetary anomalies due to central caustics in high-magnification events. Within 2 days, 
preliminary models were circulated and within 4 days, a planetary model was found that 
matched all the major lightcurve features. See Figure O 

In accord with standard microlensing practice, one person (JCY) was assigned to sys- 
tematically review all the evidence and propose a final model, which would then be vetted 
by all groups contributing data. Her report stated that the observed deviations were most 
likely due to either systematics in the data or stellar variability and so most likely implied 
that there was no planet or, in any case, that it was impossible to reliably claim a planet. 
Note that none of the evidence presented here that MOA-2010-BLG-523S is an RS CVn star 
entered into JCY's reasoning or report. 

Rather, JCY was led to question the planetary model because of three features. First, 
the model source crossing time was almost exactly half the naive time derived from inspection 
of the lightcurve. To enable this, the model has the source pass the planet-star axis almost 
exactly at right angles, so that it passes the middle (weak) cusp almost exactly at peak. 
See Figure Second, one of the three predicted features of the peak lightcurve takes place 
in a small gap in the data. Third, another feature is somewhat more pronounced in the 
model than in the data. Each of these is, by itself, quite plausible and within the range of 
microlensing experience, but together were suspicious. 

Hence, JCY sought confirmation of the planetary signal in other data sets. Unfor- 
tunately, the two other Chile observatories that might have taken such data (OGLE and 
La Silla) had very sparse coverage. She therefore investigated the CTIO ANDICAM H- 
band data, which are normally of substantially lower quality than the /-band data and so 
are usually used only for sp ecial purposes, such as comparison with high-resolution post- 



event i/-band irn aging (e.g., IJanczak et al.l 120101 ) or when /-band data are saturated (e.g. 



Dong et al.ll2009[ ). In this case, the //-band data showed a smooth peak, which is clearly in- 
consistent with the "bump" seen in /-band. See Figure |6l (Note, however, that the //-band 
peak is still asymmetrically offset by 1.5 hours compared to the time expected based on data 
in the wings.) 

It is in principle possible to have "sharper" features in /-band than //-band due to limb- 
darkening effects. This is contrary to the general expectation that microlensing is achromatic 
since in general relativity geodesies do not depend on wavelength. This exception occurs 
when the lens resolves the source because limb darkening is more severe in bluer passbands 
making the light profile more compact. Nevertheless, the amplitude of the difference seen in 
this event is much too big to be explained by this effect. The relative difference in effective 
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source sizes is 



(mm 
Mm 



1 



o.2r 



H 



1 - 0.21/ 



1 ~ 0.1(1/ -r^) <o.o2 



(14) 



whereas the difference in timescales of the observed deviations is a factor ~ 2. Hence, there 
is no plausible reason for the difference in / and H band over the peak. Moreover, JCY found 
that the other part of the "planetary signal" , the asymmetry in the lightcurve (both /-band 
and H-hand), can be fit by "xallarap" (orbital motion of the source about a companion) 
with periods of 3-15 days. 

The path to gathering the evidence summarized in Section H] was circuitous. First, in 
response to JCY's report, AU reiterated that the source (or at least some star in the aperture) 
was a variable with an 11 day period and a 1% amplitude. This variability had previously 
been ignored in the analysis due to the fact that the "planetary" deviation had a much 
shorter timescale. Then AG learned from stellar-interiors expert MHP that RS CVn stars 
were found v ery frequently at T ^ 5250 K. MHP then suggested that the Rossby number 
scalings from lHartman et al.l (120091 ) could explain the observed variability. In the meantime, 
it was found that variability at fixed period but random phase (characteristic of spots) is 
strongly favored over a strictly periodic signal. These results were consistent with an RS 
CVn star, so one would expect to see the str ong calcium H fc K emission characteristic of 
such stars in the UVES spectrum. However, iBensby et al.l ( 120111 ) had not remarked upon 
this because the blue spectral channel is rarely if ever examined for stars in this program 
because they are very heavily reddened. A check of the blue channel indeed showed strong 
H & K lines. These lines proved to be easily detectable in this case not only because they 
are intrinsically strong, but also because t he spectrum was ta ken at / ~ 13.3, which is 
substantially brighter than is typical for the iBensby et al.l (120111 ) sample (see their Fig. 1). 



In brief, the contradiction between the optical and IR lightcurves proved to be the 
crucial turning point in debunking the "planet" , even though a more detailed investigation 
of the available data provides overwhelming evidence that this was a microlensed RS CVn 
star. 

This history argues for caution in the interpretation of planetary signals, particularly 
when they are both of sm all amplitude and without the discontinuous slopes characteristic 
of caustic crossings (e.g., ICould et al.ll2006l ). One may counter in this case that RS CVn 
stars are extremely rare, but the fact remains that this "rare event" occurred within the first 
dozen or so microlensing planets. Such rare events in small samples remind us to be vigilant 
about our assumptions. 



We note that the misinterpretation of microlensed spots as planetary signals was sug- 
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gested more than a decade ago by lHeyrovsky fc Sasselovl (I2OOOI ). who specifically cautioned 
on the difficulty of distinguishing spots from planets in high-magnification events and even 
suggested intensive multi-band photometry as a means to tell the difference. As they re- 



marke d, such multi-band (optical/IR) photometry had already been advocated by lGaudi &: Gould 
(119971) as a means to bett er characterize planetary perturbations. This earlier paper (see 
also ICould &: Welch 1996) was the motivation to build the optical/IR dichroic camera 



ANDICAM f PePov et al 



20031 ). whose I/H observations of MOA-2010-BLG-523 proved 



crucial in demonstrating that the deviations were due to spots rather than a planet. There 
were several ot her early investiga tions of the interpenetration of spots and binary or planetary 
microlensing. iHan et al.l ( 120001 ) argued that spots might be easier to study in binary-lens 
than single-lens microl e nsing because the caustics were more likely to transit the source. 
And iRattenbury et al.l (120021 ) made a broader investigation of whether spots could in fact 
be mistaken for planets, arguing that this was really only possible in the rare events (such 
as MOA-2010-BLG-523) in which the lens passes very close to or over the source. However, 
to our knowledge, there are no previously published observations of microlensed spots. 



7. RS CVn Stars in the Galactic Bulge 



As we have emphasized, it will be quite rare that an RS CVn star is magnified sufficiently 
to get a high S/N spectrum of the heavily extincted C a H fc K lines . Neve rtheless, there 
are other paths toward identifying bulge RS CVn sta r s. lUdalski et al.l (120121 ) found optical 
counterparts to X-ray sources from the IJonker et al.l (120 111 ) Galactic Bulge Survey (GBS) 
catalog, including 81 spotted stars, which are very probably RS CVn stars. However, because 
the underlying X-ray catalog is confined to 1 < |6| < 2, it is likely that a large fraction of 
these are in the disk. 

There are two relatively straightforward ways to distinguish between bulge and disk 
membership. First, a subset of seven of these 81 stars are eclipsing. All but one of these 
are relatively bright 12.7 < I < 15.2, and so it should be possible to obtain spectra and 
thus measure their distances using the method of eclipsing binaries. Even the faintest of 
these, at / = 17.4, is not beyond reach. A large fraction of the remainder could be put 



on a clump-centric color magnitude diagram (INataf et al.ll201ll ). In most cases, this should 



cl early distinguish b etween bulge and foregrou nd stars. Unfo r tunat ely, the extinction map 



of iNataf et al.l (|2012| ) does not reach most of the lUdalski et al.l (|2012| ) stars because this map 
is based on OGLE-HI data, whereas the GBS survey is restricted to low-latitude fields that 
are only covered by OGLE-IV. However, it should be possible to apply the clump-centric 
method to these OGLE-IV fields as well. 
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8. Conclusions 

The evidence presented in Sections H and |5] that MOA-2011-BLG-523S is an RS CVn 
star is overwhelming. This star shows quasi-periodic variations (the form expected for spots) 
with a period of P = 10.9 days. The amplitude of variation (few percent) is consistent with 
what one would expect from its inferred Rossby number. It exhibits very strong calcium H 
& K emission, such as is seen in only 3 out of a sample of 234 local subgiants. All three 
are spectroscopic binaries, and two are known to have periods of P ~ 30 days, which given 
their radii and measured vsini implies that they are tidally spun up. MOA-2010-BLG-523S 
has high microturbulence measured compared to the 26 microlensed dwarfs and subgiants, 
particularly among stars of similar temperature. 

Unfortunately, it is not possible to say definitively whether MOA-2010-BLG-523S is in 
a binary or not. Its period relative to its radius is suggestive of being in the same class of 
tidally spun up binaries that includes the 3 calcium-active subgiants just mentioned. But its 
temperature is near the range of active subgiants found from Kepler seismology, the great 
majority of which must be single stars (or widely separated, non-interacting binaries). If a 
strong case could be made that this was not a binary, then from the lithium measurement 
(Section 14. 5p this would be evidence for an intermediate age population. But this is not 
the case. The fact that the peak is offset from the time expected from the wings by 1.5 
hours strongly suggests "xallarap" (orbital motion of the source due to a binary compan- 
ion). However, the irregular character of the lightcurve, probably due to microlensed spots, 
compromises our ability to make a rigorous microlensing fit for xallarap. 

The /-band lightcurve is well fitted by a planetary-lens model. The path to discovering 
that this is a coincidence, and that the lightcurve anomaly is due to spots was quite cir- 
cuitous, as described in Section [61 This argues for caution in the interpretation of planetary 
microlensing events in which the deviations are small and lack features that are obviously 
due to a 2-body lens. 
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Fig. 1.— Lightcurve of MOA-2010-BLG-523S from the 2001 high-cadence OGLE transit 
campaign, binned by day. There are a total of 786 observations on 32 nights, spread over a 
46 day interval. The underlying data have typical errors of 0.10 mag unbinned, hence 0.02 
mag when binned. The source shows periodic or quasi-periodic oscillations with a period of 
roughly 12 days. 
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Fig. 2. — Goodness of fit of strictly periodic (red) and quasi-periodic (black) models of 
variability of MOA-2010-BLG-523S over nine OGLE-III seasons (2001-2009). The strictly 
periodic models have 12 free parameters (period, amplitude, phase, plus zero-point offsets for 
each season), while the quasi-periodic models have 28 (additional phases and amplitudes for 
each season). The ordinate shows the difference in relative to a model with 9 parameters 
(zero-point offset at each season), [x^ = 2686.75 for 2531 dof], taking account of the different 
number of dof. Except for a spike very close to 1 day (0.9947 it 0.0005 day) the highest peak 
is at Pq = 10.914 ±0.055 day. The quasi-periodic models are clearly favored over the strictly- 
periodic ones. Other notable peaks are at the alias of the sampling frequency (0.5 day), and 
at the aliases of the main peak, P± = 1/ (I/Pq i V-^ay synod) — (0.913,1.098) day. 
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Fig. 3.— Lower panel: iBensbv et all (120 111 ) UVES spectrum of MOA-2010-BLG-523S in the 
region of the calcium H & K lines. The mean counts per 0.05 A in V and R "continuum" 
passbands are 246 and 439, respectively. Upper Panels: Zooms of the cores of the calcium 
H & K lines. These have total counts of 5596 and 5297 respectively. Hence the 5* parameter 
isS = (5596 + 5297)/[(246 + 439) /0.05 A] or 5 = 0.79 A. 
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Fig. 4. — Microturbule nce parameter £ vs. tempe rature T for 26 microlensed dwarfs and 
subgiants measured by iBensby et al.l (l2010l . 1201 ll ). "Dwarfs" {logg > 4.2), "regular sub- 
giants" (4.0 < logg < 4.2), and "large subgiants" [logg < 4.0) are shown in black, red, and 
green respectively. MOA-2010-BLG-523S is a clear outlier to the sample. 
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Fig. 5.— Planetary model of MOA-2010-BLG-523 (black) fit to /-band data points from 
several observatories as indicated in legend. i/-band data are not shown. 
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Fig. 6. — Residuals to a point-lens (finite source) fit that uses only iJ-band data over the 
peak, for CTIO iJ-band (crosses) and /-band (circles). Points are binned in 10- minute 
intervals. Error bars (not shown) are slightly smaller than points. In contrast to the /-band 
data, the //-band data show no convincing case for strong deviations. 



